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Abstract. The first part of this paper defines recursive interac- 
tions by means of logistic functions and derives a general result on 
the way interacting systems evolve in attractors. It also defines 
the notion of coevolution trajectory and presents a new family of 
attractors: orbital attractors (including single, irregular, folded, 
complex and discontinuous orbits). The second part summarizes 
the results of a first experimental analysis of recursive interactions 
in both binary and multiple interactions. Among other results, 
this analysis reveals that interacting systems may easily evolve 
from chaos to order. 

1. Introduction 

Mathematical continuity makes no sense in biology. The biosphere is 
radically discontinuous, i.e. discrete; it is the paradigm of discreteness. 
Each living being is a discrete unity, unique and unrepeatable, that 
emerges and self maintains at the expense of a discrete net of metabolic 
reactions governed by a discrete net of discrete information units. Each 
living being forms part of a discrete group of discrete individuals which 
in turns form part of a discrete net of ecological interactions. But living 
organisms are not the only discrete objects in the universe, matter and 
energy are also discrete. In consequence, all natural processes involving 
matter or energy interchanges have to be of a discrete nature. Even 
space and time could be discrete as has been repeatedly suggested from 
different areas of physics ([6], [7], [9], [19], [13], [10], [24], [4], [24], [2], 
[21], [11], [5], [3], [22], [1] , [12], [26], [23], [17], [3]). In these conditions, 
and being the best model of any object the very object itself, what 
should be put into question is not the role of discrete models [8] but the 
role of mathematical continuity in an essentially discontinuous world. 

In addition of being discrete and interactive, living organisms are also 
recursive. Each generation determines the next one through a repro- 
ductive process that is not exact, giving therefore the appropriate op- 
portunity to evolution. It is then clair the reason for which biology has 
always been interested in discrete and recursive models. Among them 
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the logistic function, perhaps the most simple and productive model of 
some biological significance. From the pioneering works of Stanislaw 
Ulam, Paul Stein [25], Nicholas C. Metropolis [20] and Robert May [18] 
to nowadays, iterative calculus and the logistic function have been an 
inexhaustible source of surprising results, from deterministic chaos to 
periodic attractors; results that, on the other hand, were immediately 
generalized by the so called principle of universality. As is well known, 
the control parameter A of a logistic function: 



determines the system evolution, particularly the type of attractor it 
finally falls in. Now then, biological systems do not evolve separately 
but immerse in a complex interaction network. For this reason biol- 
ogy has also been interested in modeling interactions. Most of those 
models, as the classical Lotka-Volterra, are defined in terms of differ- 
ential (or finite differences) equations. But as far as I know, biological 
interactions have never been modeled by logistic functions. As we will 
immediately see, biological (and non biological) interactions can easily 
be modeled by means of this type of functions. And the results can- 
not be more interesting from both the mathematical and the biological 
point of view. In fact, on the one hand recursive logistic interactions 
are a new source of new mathematical objects as coevolution trajecto- 
ries or orbital attractors. On the other, they provide us a new way 
of examining coevolution processes and derive some significant results 
on the way systems coevolve. For instance, that interacting systems 
may coevolve from chaotic regimes to stable states defined by peri- 
odic attractors of low period; or that complex attractors, as fractal 
or chaotic, may suddenly evolve in single periodic attractors; or that 
systems to the brink of extinction may be completely recovered and 
stabilized thanks to its recursive interactions with other systems; or, in 
the case of multiple interactions, that system stability grows with the 
number of interactions. 

The discussion that follows begins by defining the appropriate theo- 
retical notions, from which a general result related to the way interact- 
ing systems evolve in attractors is formally derived. In the following 
two sections I define some new mathematical objects as the coevolution 
trajectory of two interacting systems and a new family of attractors. 
Finally, I resume the results of a first experimental work on recursive 
logistic interactions, including a short presentation of multiple interac- 
tions. Evidently, the field to explore is immense and what follows is 
but a simple introductory note. All numerical and graphical data have 
been obtained with the aid of InterCalculus [16], a computer applica- 
tion developed to analyze recursive interactions by means of logistic 
functions. 
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2. Recursive logistic interactions: Definitions 

Let us consider two systems X and Y modeled by two logistic real 
functions: 



An elemental way to make both functions mutually dependent consists 
in defining the control parameter of each function in terms of the other 
function variable: 



If gxiVn), the interaction function of system X with system Y, is de- 
creasing the interaction of system X with system Y will be said neg- 
ative; if it increasing the interaction will be said positive. The same 
applies to gy{xn), the interaction function of system Y with system 
X. There are many ways to define the interaction functions and gy. 
Here we will use the following: 

gxiVn) = 4(1 - S^yn) 

for negative interactions, and: 

for positive ones. In both cases, represents the sensitivity of system 
X to system Y , i.e. a measure of the effects of system Y on system 
X. It will take its values within the real half closed interval (0,1]. 
Similarly, for the interaction of system Y with system X we will have: 



It is now immediate to define the following three types of recursive 
logistic interactions (recursive interactions from now on): 

(1) Negative-Negative: 




■^n+l gx{yn)-^n{^^ •^n) 
Vn+l = gy{Xn)yn{l " Vn) 



gy{Xn) = 4(1 - SyXn) 
gyi^Xn) 4:SyXn 



Vn+l 



4(1 - S^yn)Xn{l - Xn) 
4(1 - SyXn)yn{l - Vn) 



(2) 
(3) 



(2) Positive-Negative 



Vn+l 



Xn+l 



^SxynXniX Xri) 

4(1 - SyXn)yn{^ - Vn) 



(4) 
(5) 



(3) Positive-Positive: 



yn+l = 4:SyXnyn{l - yn) 




(6) 
(7) 
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Here we will exclusively deal with these recursive interactions. From 
now on negative-negative interactions will be referred to as NN inter- 
actions; positive-negative as PN; and positive-positive as PP. 

The above definitions can immediately be generalized for any num- 
ber of systems. In effect, consider m systems (^j)i=i,2,...m being each 
modeled by a logistic function of a real variable Xi. Assume that all 
systems interact which each other. Their evolution can be expressed 
as: 

^ m 

■^j, n+1 ^ fl'i^i •^i,n) (8) 

where Xi^n represent the n-th generation of system Xi] gxi{xj) is the 
interaction function of system Xi with system Xj, a function that will 
be given by iSijXj^n if that interaction is positive, or 4(1 — SijXj^n) 
if it is negative, being Sij the sensitivity of system Xi to system Xj. 
Equations (8) can be easily computerized so that it is possible to make 
both numerical and graphical analysis of the evolution of any number 
of interacting systems. 

3. Attractors theorem 

We will prove now a general and basic result that states the simultane- 
ity and equivalence of the attractors the interacting systems evolve in. 

Theorem 1. Let X andY be two systems that interact with each other. 
A necessary and sufficient condition for system X to fall in a period-k 
attractor is that system Y also falls in a period-k attractor. 

Proof. Assume that system X falls in a period-k attractor. From a 
certain integer i > it will hold: 

Xn = Xn+k, Vn > i (9) 
On the other hand we will have 

Xn = 5'(y„_i)x„_i(l - Xn~l) 
Xn+k = S(jjn+k~l)Xn+k-l{^ ~ Xn+k-l) 

where S(yn-i) is 4:{l — Sxyn-i) if the interaction of X with Y is negative, 
or ASxHn-i if it is positive. The same applies to S{yn+k~i)- According 
to (9) we can write: 

S{yn-l)Xn-l{l - = S{yn+k-l)Xn+k-l{'^ " Xn+k-l) 

and taking into account that Xn-i = Xn+k-i we have: 

S{yn~i) = S{yn+k^i) 

That is to say: 



4(1 - S^yn-i) = 4(1 - S^yk+n-i) 
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if the interaction of X with Y is negative; or: 

if it is positive. In both cases, we immediately get: 

Un-i = Vn+k-i, Vn > i 

which imphes that system Y has also fallen in a period-k attractor. 
This proves that if system X evolves in a period-k attractor, so evolves 
system Y. A similar argument proves the complementary result. 

□ 

The above theorem applies only to binary interactions. In the case 
of multiple interactions, in fact, systems do not reach simultaneously 
their final attractor; although, in most cases, once a first system falls in 
its final attractor the remainder ones also fall in their respective final 
attractors after a few number of iterations. I think attractors theorem 
is not the only conclusion that can be formally derived for binary inter- 
actions. Experimental research strongly suggests that, among others, 
the following results could also be formally proved: 

• In NN interactions, the less sensitive system always goes ahead^ 
of the most sensitive one. 

• In PN interactions, the system suffering the negative interaction 
always goes ahead of the other. 

• Coevolution trajectories in NN interactions are symmetrical 
with respect the bisectors of (0, 1) x (0, 1). 

4. Coevolution trajectories 

Recursive interaction invariably leads to a final attractor that may be 
periodic, chaotic, fractal or of some new types, as we will immediately 
see. The set of points each interacting systems traverses towards its 
final attractor defines its coevolution trajectory (CT). It usually con- 
sists of two or more well defined lines, branches from now on, usually 
with diagonal symmetry although they can also exhibit radial or spiral 
symmetry. Each interacting system has its own CT, being both CT 
very similar. 

Systems evolve along the branches of their CT by cyclically step- 
ping from a branch to other, always in the same order, so that they 
exhibit a remarkable periodic behaviour. Due to the discrete nature of 
iterations, CT are not continues but discrete lines, although the den- 
sity of points can be extremely high and we have to magnify them 
thousands of millions times to make discontinuity visible. Naturally, 
in these superdense regions the progress towards the final attractor is 
also extremely slow. Coevolution trajectories are then sets of points 



System X goes ahead of system Y if x ^ y. 
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-usually lines- in the interval [0, 1] x [0, 1] of along which interact- 
ing systems evolve towards their respective final attractors. As far as 
they have been examined, the coevolution trajectories of both systems 
are geometrically similar, with the same number of branches and the 
same type of symmetry. The branches of a CT may converge or diverge 
from the central point (0.5,0.5), being the changes simultaneous in all 
branches of both CT. 

There also exists CT trough which systems evolve from chaos to non 
chaotic attractor (periodic, fractal or orbital) or even from a chaotic 
regime to other different chaotic regime. While traversing their re- 
spective CT, one of the systems goes always ahead of the other (the 
less sensitive in the case of NN interactions, and the one suffering the 
negative interaction in PN interactions). With the appropriate initial 
seeds it is possible, however, to start from an inverted position. In 
these cases an inversion of the CT occurs thanks to which the systems 
recover its normal relative position. To go ahead of the other seems to 
be a general law that operates even in chaotic regimes. 

5. Orbital attractors 

As in the case of single logistic recursion, recursive interactions also end 
up by reaching a final attractor. Apart from the well known periodic, 
fractal and chaotic attractors, in the case of recursive interactions we 
can also observe at least a new family of attractors that we will term 
orbital attractors. An orbital attractor is typically composed of one or 
more closed lines (orbits) which are usually eccentric. Each orbit its 
initially a discontinuous set of segments which progressively extend in 
the same direction so that finally they overlap and close the line. In 
some cases the orbits are single lines while in other they have a com- 
plex internal structure. Systems trapped in an orbital attractor behave 
as if they were orbiting, although their true behaviour is a little more 
complex. In fact, they jump cyclically from an orbit to other always 
in the same order and in such a way that the successive jumps on the 
same orbit take place on different segments, always in the same order. 
All segments grow in the same direction so that after a few interactions 
each segment reaches the next one, although, being of a discontinuous 
nature, their respective points do not coincide, each segment continues 
indefinitely by occupying the empty gaps along the orbit. After a con- 
siderable number of iterations, each segment reaches itself and progress 
trough its own gaps. This way of progressing suggest the possibility 
that after a huge number of iterations all orbital attractors finally end 
in a periodic attractor. 

The orbits of an orbital attractor may be at least of one of the 
following types: 

(1) Single 

(2) Irregular 
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(3) Folded 

(4) Complex 

(5) Discontinuous 

The lines of a single orbit do not exhibit internal structure, they always 
appear as single lines even if we magnify them by several billions times. 
The same apply to irregular and discontinuous orbits although in these 
cases the lines exhibit more or less irregular forms, as complex indenta- 
tions. Folded orbit exhibit a complex non fractal internal structure of 
folded lines in the regions of maximum curvature, this structure seems 
disappear in the less curvature regions. This is also the case of com- 
plex orbits, although in these orbits the internal structure is extremely 
complex, perhaps of a fractal nature. In the case of discontinuous or- 
bits, the segments do not extend and maintain the empty gaps between 
them. After a certain number of iterations these attractors degenerates 
in a single period attractor. 

A remarkable characteristic of attractors in recursive logistic interac- 
tions is that chaotic, fractal and orbital attractors may suddenly evolve 
in a periodic attractor whose points presumably belong to the original 
attractor. It is also remarkable the existence (at least in NN interac- 
tions) of periodic attractors whose attraction decreases exponentially 
as systems approach them, so that it takes thousands of millions of iter- 
ations to progress one decimal cipher towards the attractor final value. 
They could be termed asymptotic attractors. It is possible that no fi- 
nite number of iterations suffices for the system to attaint the attractor 
(in the same sense that the limit of a sequence cannot be reached by 
the successive terms of the sequence). 

6. Negative-negative interactions 

In accordance with the above definition (2)-(3), in NN interactions the 
growth of a systems is always to the detriment of the other. Thus, 
NN interactions model the coevolution of two systems that compete 
with each other. The results will depend on the sensitivities and on 
the initial seeds of both systems, although the sensitivity dependence 
is stronger. The coevolution of interacting systems is therefor con- 
trolled by four real number in the interval (0, 1], which means a huge 
number of possibilities to examine, each representing a possible coevo- 
lution history. Although only an insignificant number of cases has been 
examined, the following conclusions could be of general application: 

(1) The branches of the coevolution trajectories are symmetrical 
with respect to the bisector, Xn+i = x„ in the case of system 
X, or yn+i = Un in the case of system Y). 

(2) All branches of a CT are geometrically similar, with same length 
and the same point density. 
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(3) Systems traverse the branches of their respective CT by step- 
ping from one branch to other, always in the same order so that 
they exhibit a remarkable periodic behaviour while evolving to- 
wards the final attractor. 

(4) The oscillations in NN interactions are synchronic, i.e. low and 
high values are simultaneously reached by both systems at the 
same successive iterations. 

(5) There is a high degree of correlation between the successive 
values reached by both systems. The correlation coefficient is 
in most cases greater than 0.9, even if both systems have evolved 
in chaotic attractors. 

(6) The system of less sensitivity go always ahead of the other. 
That is to say, the values successively reached by the less sensi- 
tive system are always greater than those reached by the other. 
This result coincides with St Matthew Theorem, a formal con- 
clusion derived from independent thermodynamical considera- 
tions (internal entropy production) [15]. 

(7) Systems evolves in an attractor that may be periodic, chaotic, 
fractal, orbital or asymptotic. Low period attractors are per- 
haps the most frequent in NN interactions. 

(8) Both systems evolves always in the same type of attractor. 

(9) Each point of a periodic attractor has its own independent 
branch in the CT. Although periodic attractors can also result 
from the degeneration of a complex attractor (chaotic, fractal 
or orbital). 

(10) When systems fall in orbital attractors, the orbit of the less sen- 
sitive system always raises over the orbit of the more sensitive 
one. It is also of less size, which means that the less sensitive 
system is more stable. 

(11) Systems can evolve from a chaotic regime to a single periodic 
attractor. 

(12) Systems can evolve from order to chaos after traversing their 
respective regular and symmetrical CT. Or in other words, after 
a long stage of periodic behaviour the interacting systems can 
evolve to a chaotic regime. 

(13) Systems can evolve from an initial chaotic regime to other dif- 
ferent final chaotic regime. Between both chaotic regimes, sys- 
tems exhibit a regular periodic behaviour while traversing their 
respective CT. 

(14) Chaotic attractors extend on a broad region and exhibit a typ- 
ical internal structure of parabolic gaps. 

(15) In certain cases, both systems have similar CT although they 
are traversed in opposite senses (trajectory inversions). 

(16) The extinction of a system (to reach the value of 0) is extremely 
rare. 
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7. Positive-negative interactions 

According to (2)-(3) the effects of PN interactions are different in both 
interacting system: one of them benefits from the growth of the other 
while this other is negatively affected by the growth of the first. This 
type of interaction, therefore, models prey-predator interactions. As 
in the case of NN interactions, the coevolution of both systems also 
depends on four real variables within the real interval (0, 1): the sensi- 
tivities and the initial seeds of both systems. Consequently there is also 
here a huge number of possibilities to explore. The behaviour diversity 
is now even greater than in the NN case. Some of the most relevant 
characteristics of this type of recursive interaction are the following: 

(1) The coevolution trajectories of both systems are similar, al- 
though the length of the branches may be different. 

(2) The coevolution trajectories can be convergent or divergent, or 
may suddenly change from convergent to divergent (or vicev- 
ersa). As in the NN case, they are symmetrical with respect to 
the bisector Xn+i = Xn in the case of system X, or = i/„ in 
the case of system Y. 

(3) The system suffering the positive interaction (predator system) 
is more affected by the interaction and may evolve to extinc- 
tion. In these cases the other system (prey system) remains in 
a chaotic attractor. 

(4) It has not been observed that both systems evolves in a chaotic 
attractor. 

(5) Both systems may evolve from a chaotic regime to any other 
non chaotic attractor. 

(6) Being on the brink of extinction, the predator systems may 
suddenly recover and evolve to an stable non chaotic attractor. 

(7) For high values of sensitivities, the branches of both CT are 
either radial or spiral. 

(8) Spiral trajectories may be more or less dense, and the number 
of their corresponding branches is variable. 

(9) The center of spiral trajectories is a single period-1 attractor. 

(10) As in the cases of natural prey-predator interactions, systems 
oscillate asynchronously, so that when a system reaches a high 
value the other reaches a low one, and viceversa. 

(11) Contrarily to the NN case, trajectory inversions have not been 
observed. 

(12) Single attractors of period-1 are frequent. They may be reached 
through one or more (convergent) branches. 

(13) Some orbital attractors exhibit very complicated forms in their 
regions of maximum curvature (folded orbits and complex or- 
bits). 
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(14) Orbital attractors of the systems suffering the negative interac- 
tion always raises over the orbital attractor of the system suf- 
fering the positive interaction. They are also of less size than 
it. 

8. Positive-positive interactions 

PP interactions model pure cooperation. There are sufficient exper- 
imental reasons to conclude that in all cases both systems evolve to 
extinction. The only exception is the elemental case Sx — Sy — Ij 
= Ho = 0.5, which evolve to the single attractor 0.5. Even in the 
case of multiple interactions systems evolve to extinction if the number 
of PP interactions is sufficiently high. From thermodynamic analysis 
we know, on the other hand, that this type of interaction does not 
have asymmetrical coevolution trajectories of minimum entropy pro- 
duction. The only trajectory of minimum entropy in this case is the 
bisector y = x. 

9. Multiple interactions 

Nature is not a set of isolated couples of systems that interact with 
each other. It is rather composed of an immense and complex inter- 
action network in which participate millions of systems that interact 
with many other different systems. It makes sense therefore to con- 
sider the possibility of analyzing more complicated situations than the 
simple binary cases we have just examined. As we have seen, binary 
logistic interactions can be easily generalized to any number of inter- 
acting systems (8). In these conditions, attractors theorem no longer 
holds and the theoretical study is much more complicated than in the 
binary case. Fortunately, we can make use of computer applications to 
analyze the evolution of complex interaction networks. 

With the aid of the application above mentioned, some experimental 
work have been performed involving up to 1000 interacting systems, be- 
ing each system capable of interacting (positively and negatively) with 
a variable number of other systems, including the case of all the others. 
For the sake of clarity, let us term stable to those systems that have 
evolved in a very low period attractor (usually a single attractor); and 
network stabilizing capacity to the number 1/r where r is the number 
of iterations which are necessary for all systems to become stable. In 
these conditions, the main and more remarkable conclusion we immedi- 
ately get from the above experimental analysis is the great stabilizing 
capacity of complex interaction networks. By way of example, in a 
set of 1000 systems in which each system interacts positively with 100 
systems and negatively with other 400, all systems become stable in 
less than 40 recursive logistic interactions. Next, I resume the most 
interesting conclusions of the performed experiments: 
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(1) The stabilizing force of the interaction network grows as the 
number of interacting systems grows; and for a definite number 
of systems, as the number of interactions per system grows, and 
as the number of PN interactions grows. 

(2) Collective behaviour has been observed in some interactions. 
In these cases systems oscillate synchronically according to a 
periodic pattern. 

(3) The standard deviation of the network stabilizing capacity de- 
crease as the number of systems increases and as the number of 
interactions per system increases. 

(4) The number of stable systems increases quickly once the first 
system becomes stable. 

(5) The final average Xj of all systems depends on the particular 
experiment, but is always around a value of 0.6. This number 
decreases towards 0,499 as the number of positive interactions 
increases. 

(6) The standard deviation of Xj is always very low, usually less 
than 0.00015. 

(7) Complex networks of PN interactions are extremely stable. The 
average value of systems in this case is 0,499, with a standard 
deviation less than 0.00007. 
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The following figures are examples of coevolution trajectories and at- 
tractors in both NN and PN interactions. Points {xn,Xn+i) of system 
X and points i/n+i) of system Y are plotted in red and blue respec- 
tively; attractors are plotted in black (for more information and figures 
download [14]) 




Figure 1: Spiral CT in a case of PN interaction {S^ = 0, 999999; Sy = 
0,99; xo = 0,9; yo= 0,9; zoom: 28419). 




Figure 2: Left: From Chaos to order in a NN interaction [S^ = 0,9998; 
Sy = 0,999; Xo = 0,001; yo = 0,9). Initially both systems evolve in a 
chaotic regime -parabolic regions- and then, after some thousands inter- 
actions, they evolve through their respective symmetrical coevolution 
trajectories towards a final 2-period attractor. 
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Figure 3: Orbital attractor of single orbit in he interaction defined by 
= 0,988; Sy = 0,3; Xq = 0,9; yo = 0,9; zoom: 4. 
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Figure 4: Orbital attractor of irregular orbit. NN interaction = 
0,335; 5*^, = 0,333; Xq = 0,3; yo = 0,2; zoom: 3. 
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Figure 6: Orbital attractor of folded orbit corresponding to the PN 
interaction Sx = 0,999; Sy = 0,232; xq = 0,9; yo = 0,9; zoom: 3 
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Figure 7: Complex region of an orbital attractor of folded orbit corre- 
sponding to the PN interaction = 0,999; Sy = 0,232; Xq = 0,9; i/q = 
0,9; zoom: 97. 
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Figure 8: Diagonal symmetry towards a final period- 16 attractor corre- 
sponding to the NN interaction defined by = 0.99988; Sy = 0.9976; 
Xo = 0.765; yo = 0.234; zoom: 2. 



16 



BEYOND CHAOS 




Figure 10: Radial CT towards an orbital attractor of = 0.999; Sy = 
0.393999; xq = 0.45; yo = 0.67; zoom: 7 
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Figure 11: Radial CT towards a period- 1 attractor of the PN inter- 
action defined by = 0.999; Sy = 0.399; Xq = 0.8; yo = 0.8. zoom: 
115. 
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